Power-Law tailed statistical distributions and Lorentz transformations 
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The present Letter, deals with the statistical theory [Phys. Rev. E 66, 056125 (2002) and Phys. 
Rev E 72, 036108 (2005)], which predicts the probability distribution p(E) oc exp^(— /), where, 
I oc PE — P/i, is the collision invariant, and exp^(a;) = (Vl + k.'^x^ + kx)^^'^, with < 1. This, 
experimentally observed distribution, at low energies behaves as the Maxwell-Boltzmann exponential 
distribution, while at high energies presents power law tails. Here we show that the function exp^ (x) 
and its inverse lnK(a;), can be obtained within the one-particle relativistic dynamics, in a very simple 
and transparent way, without invoking any extra principle or assumption, starting directly from the 
Lorentz transformations. The achievements support the idea that the power law tailed distributions 
are enforced by the Lorentz relativistic microscopic dynamics, like in the case of the exponential 
distribution which follows from the Newton classical microscopic dynamics. 

PACS numbers: PACS number(s): 51.10.+y, 05.20.-y, 52.27.Ny 



I. INTRODUCTION 

The origin of the power-law tails in the particle pop- 
ulation, observed in a wide variety of systems from high 
energy physics to condensed matter physics, is one of the 
most pressing outstanding issues in statistical physics. 
For instance the cosmic ray spectrum, which extends 
over 13 decades in energy, from 10* eV to 10^° eV, and 
spans 33 decades in particle flux, from 10~^^ to 10^ 
units, is not exponential at relativistic energies [ll-Q- 
More precisely, the particle spectrum fi oc xiP^i ~ /^m) 
obeys the Boltzmann law of classical statistical mechan- 
ics i.e. x(a;) ^ exp(— x) for low energies, while for high 



energies this spectrum presents power law fat tails i.e 



X{x) 



■ -t-oo 



the spectral index being close to 2.7-3.1. 



Also theoretical problems present difficulties when 
treated within the classical statistical mechanics. For in- 
stance the Bekenstein-Hawking area law problem, related 
to the black hole entropy, asserts that the entropy of a 
black hole scales as the area of the event horizon [4] . It is 
well known that the ordinary Boltzmann statistical me- 
chanics and thermodynamics fail to justify this law, and 
conversely predict that the entropy scales as the volume 
of the spatial region delimited by the event horizon. 

In experimental physics, the cosmic rays problem was 
approached for the first time in 1968, by using a differ- 
ent distribution from the Boltzmann one. In his proposal 
Vasyliunas heuristically identified the function x{^) with 
the Student distribution function which presents power- 
law tails 0]. In the last 40 years a vast amount of lit- 
erature has been produced, based on the Vasyliunas dis- 
tribution, including both experimental and theoretical 
works regarding statistical systems violating the Boltz- 
mann statistics. 



On the other hand, in the last few decades, the power- 
tailed statistical distributions have been observed in a 
variety of physical, natural or artificial systems. In Q 
one can find updated discussion of power-law tailed dis- 
tributions and an extensive list of systems where that 
distributions have been empirically observed. 

In several papers specific physical mechanisms have 
been explored in order to furnish theoretical support 
to the experimentally observed power-law-tailed distri- 
bution functions. However there is currently an intense 
debate regarding the theoretical foundations of the non- 
Boltzmannian distributions, especially in physics. In the 
last years, after noting that the power-law tails are placed 
in the high energy region, and then regards relativistic 
particles, the question has been posed whether the so- 
lution of the problem, i.e. the theoretic determination 
of the function x{^) ^-nd consequently of the related dis- 
tribution and entropy, can be explained by invoking the 
basic principles of special relativity. 

The present Letter deals with the statistical theory, 
predicting for the function x.{x) the very simple form i.e. 
X{x) = exp„(x) with 



exp^(x) = + + KxY^'^ . 



(1.1) 



The parameter < 1 is the reciprocal of light speed, in 
a dimensionless form, while is the particle rest en- 
ergy. In the classical limit k — 0, exp^(a;) reduces to the 
ordinary exponential function, defining the Boltzmann 
factor of classical statistical mechanics 0-01 • 

Starting from the inverse function of exp^(a;), namely 
ln„(x) given by 



ln„(x) 



2k 



the entropy of the system assumes the form 
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(1.2) 



(1.3) 
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The above entropy, after maximization under the con- 
straints expressing the conservation of the distribution 
norm and of the system energy, yields the following ex- 
pression for the distribution function 



1 



/ = -exp,(-7/? [E-^i]) 



(1.4) 



where E is the particle energy, (3 and are the Lagrange 
multipliers, while 



1 



e = exp^(7). 



(1.5) 



are two constants, depending exclusively on the value of 
the parameter k. 

In the last few years various authors have considered 
the foundations of this statistical theory, e.g. the EE- 
theorem and the molecular chaos hypothesis Q , the 
thermodynamic stability [H, [Tl|, the Lesche stability 
p^ - fTsj , the Legendre structure of the ensued thermody- 
namics [H [l3| , the geometrical structu re jisl , and other 
mathematical properties of the theory [1^, etc. On the 
other hand specific applications to physical systems have 
been considered, e.g. the cosmic rays [ij, relativistic [23| 
and classical t21i] plasmas in presence of external elec- 
tromagnetic fields, the relaxation in relativistic plasmas 
under wave-particle interactions [23l | , anomalous dif- 
fusion [2J, 25], kinetics of interacting atoms and photons 
[25} , particle kinetics in the presence of temperature gra- 
dients 1271 , particle systems in external conservative force 
fields [2I], stellar distributions in astrophysics [29l - l3]| . 
quark-gluon plasma formation [s^, quantum hadrody- 
namics models [33,], the fracture propagation [sj], etc. 
Other applications regard dynamical systems at the edge 



3£ I , the random ma- 
42| . the game theory 



of chaos [35' -'sTi , fractal systems [H, 
trix theory [40, 41j , the error theory 

[isj . the information theory Q, etc. Also applications to 
economic systems have been considered e.g. to study the 
personal income distribution [45-48], to model determin- 
istic heterogeneity in tastes and product differentiation 
gi, iOl etc. 

In the present contribution we reconsider critically the 
physical origin of the above mentioned statistical theory. 
Our main goal is to show that the function exp^(a;), and 
its inverse \n^{x), can be obtained within the one-particle 
relativistic dynamics, in a very simple and transparent 
way, without invoking any extra principle or assumption, 
starting directly from the Lorentz transformations. 



II. LORENTZ TRANSFORMATIONS AND 
K-EXPONENTIAL 

In the present section we will show that the function 
exp^(x) emerge within the special relativity as the rela- 
tivistic generalization of the ordinary exponential of clas- 
sical physics. Special relativity is a theory founded on 
the Lorentz transformations, so the question whether it 



is possible to obtain exp^(x) starting directly from the 
Lorentz transformations and how, arises naturally. 

Let us consider in the one-dimension frame S two iden- 
tical particles A and B, of rest mass m. We suppose that 
the two particles move with velocity Va and Vb respec- 
tively. The momenta of the two particles are indicated 
with Pa = p{va) and Pb = p{vb), while their energies 
are indicated with Ea = E (va) and Eb = E (vb) respec- 
tively. 

We consider now the same particles in the rest frame S' 
of particle B. In this new frame the velocity, momentum 
and energy of the particle B are v'^ — O^p'^ =0 and E'^ = 
mc^ respectively. In S' the velocity of particle A is given 
by the formula — {va — Wb)/(1 — VaVb/c^) defining 
the relativistic velocity composition law, which follows 
directly from the kinematic Lorentz transformations. In 
the same frame S' the momentum and energy of particle 
A are given by the dynamic Lorentz transformations 



p'a = 1{Vb)Pa ■ 

E'^ = ^{yB)EA 



C ^VbJ{Vb)Ea 
- Vb1{Vb)Pa , 



(2.1) 
(2.2) 



i{vb) = (1 — v%/c^)~'^/'^ being the Lorentz factor. 

We note that the expression of the momentum Pb — 
mvBj{vB) and of the energy Eb — rac^^{vB) of the 
particle B appear explicitly in the latter transforma- 
tions. After introducing in place of the dimensional 
variables {v,p,E) the dimensionless variables {u,q,£) 
through v/u = p/mq — \J EjraE = \k\c = < c, the 
Lorentz transformations become 



9a = K^^A^B - K^fe^A 

E — n EaEb ^ qAQs , 



(2.3) 
(2.4) 



while the classical limit c — > 00 is replaced by the limit 
K — > 0. From the Lorentz invariance it follows easily the 
energy- momentum dispersion relation k*E'^ — K^q^ = 1. 

From its definition it follows that k is the Einstein j3 
factor related to the velocity and the condition < c 
implies that —1 < k < -1-1. For a particle at rest it 
results £;(0) = mc^ and then 5(0) = I/k^. Then I/k^ 
represents the dimensionless rest energy of the particle. 
Alternatively 1/ n is the refractive index of a medium in 
which the light speed is f * . 

It is remarkable that the dynamic Lorentz transfor- 
mations form a system of two coupled linear equations, 
having the important feature that the contributions of 
the two particles A and _B, appearing in the right-hand 
side of the equations, are not factorized. Spontaneously 
the question emerges at this point, whether new variables 
exist, able to factorize the contribution of the two par- 
ticles A and B in the right-hand side of the dynamic 
Lorentz transformations. For the kinematic Lorentz 
transformations this problem was first posed and solved 
by Minkowski, by introducing light cone variables. Af- 
ter observing that both kinematic and dynamic Lorentz 
transformations have the same structure, it is easy to ver- 
ify that the two dynamic light cone variables are given 
by K^E zL Kq > 0. Clearly any power of these variables 
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permits the factorization of the right-hand side of the 
Lorentz transformations and therefore we can write the 
more general form of the dynamic Minkowski light cone 
variables as follows 

X±^{n''£±KqY . (2.5) 

Lorentz transformations (|2.3p and (|2.4p . after using the 
variables x±, assume the following required form 

X±A = Xts X±a , (2.6) 
while the condition expressing the Lorentz invariance of 
the system simplifies to X-X+ = 1- 

By taking into account the energy-momentum disper- 
sion relation n^E"^ — n^q^ ~ 1, the light cone variables 
X± can be viewed as functions of the momentum q i.e. 
X± — xi^l)- III order to completely define the func- 
tion xil) need to fix the exponent a appearing in 
(|2.5|) . Only if we pose a = the Lorentz transforma- 
tions can be reduced to the Galilei transformations, in 
the classical k — > limit. Indeed, for a = it results 
limK-fO xil) — cxp((7) and the Lorentz transformations 
reduce to 



exp(±g^) = exp(=FgB)exp(±g^) , (2.7) 

so that the Galilei transformation for the momenta, q'^ = 
qA — qB follows immediately. 

After posing a = 1/k, and taking into account the 
dispersion relation £ — \/l + K^q'^/n^, Eq. (j2.5l) univo- 
cally defines the dynamic Minkowski light cone vari- 
ables xil): which interestingly, coincides with the k- 
exponential given by (jl.ip . i.e. 



X{q) = exp^(g) 



(2., 



Finally the Lorentz transformations (|2.6|) assume the 
required, factorized simple form 



exp„(±g^) = exp^(=FgB)exp„(±g^ 



(2.9) 



It is remarkable that the function exp^(a;), defining 
the particle distribution function in statistical physics, 
emerges also in one-particle special relativity as the light 
cone variable, able to factorize the Lorentz transforma- 
tions, and represents the relativistic generalization of the 
ordinary exponential. 

The Galilei relativity principle, holding both in classi- 
cal physics and in special relativity, imposes the equiva- 
lence of all the inertial frames. According to this prin- 
ciple, the inverse Lorentz transformations must have the 
same structure of the direct transformations (|2.9p except 
for the substitutions q'^ •(-> and qs — > —qs- This 
requirement implies for the generalized exponential the 
property 

exp^(a;)exp„(-a::) = 1 , (2.10) 
which in terms of the generalized logarithm becomes 



(2.11) 



The latter relationships imposed directly by the Galilei 
relativity principle, represents two important properties 
of the functions exp^(a;) and \iifi{x) respectively. 



III. LORENTZ TRANSFORMATIONS AND 
K-LOGARITHM 

In the present section we will show that the function 
Iuk (x) , can be obtained by employing the Lorentz trans- 
formations only. 

Let us consider the Lorentz transformations in the case 
where q^ = q^ q^ = q -\- dq and q'^ = d^q. Eq. p.9p 
becomes 



exp^(d«(7) ^exp^{-q)e^p^{q + dq) 



(3.1) 



By performing a Taylor expansion, up to the first order, 
in the latter equation and after taking into account ()2.10p 
and the normalization condition exp'^(O) = 1 we obtain 



df,q 



= cxp«,(g) 



(3.2) 



Preliminary, we observe that the composition law of 
the relativistic momenta can be established starting from 
the first equation of the Lorentz transformations (|2.3p . 
after elimination of the energy by taking into account 
the dispersion relation £ = \/l + n'^q^ /k^, obtaining 

(3.3) 



The expression of the K-differential d^q, can be ob- 
tained starting directly from the latter composition law, 
by posing qg = q, qA = q + dq and q'^ ^ d^q and after 
performing a Taylor expansion i.e. 



df,q = 



dq 
7k('7) 



(3.4) 



7k((?) = ■\/l + K^g^ being the Lorentz factor. After tak- 
ing into account (|3.4p . Eq. p.2p becomes 



(3.5) 



, , d exp g 
7k (g) — ^ — = exp^g 



The differential equation for the function luK^f) can be 
obtained by performing in Eq. (j3.5l) the change of vari- 
able q = ln„(/). 



d In^/ _ 7k (Ins/) 



df 



f 



(3.6) 



At this point after posing exp^(g^) = /, exp^(g^) — 
/', and exp^(— g^) = e, the Lorentz transformations in 
the factorized form (|2.9p . become /' = ef. On the other 
hand the first of the Lorentz transformations (|2.3p yields 

ln,(e/) =74ln,e) In, / + 7, ( In, /) In, e , (3.7) 

so that we obtain that 7, ( In, /) is given by 

7K(ln,/)=i^-^^4^1n,/. (3.8) 
In, e In, e 

Finally, by accounting p.Sp in Eq. p.6p we obtain 
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IV. CONCLUSIONS 



d . ln.(e/) 
— / ln„/ = — 



After choosing w» 



(3.9) 



suits = — i/Vi — K 

The latter relationship simplifies Eq 



—v^ , which implies — 
^ and therefore In^ e = 1 
3:9 



- 1 , it re- 



df In^ e 



1-k2. 
as follows 

(3.10) 



Eq. p.lOp . already known in the literature, has been 
introduced heuristically in [l[ . Its solution, with the con- 
dition (|2.11l) . is unique and defines the generalized log- 
arithm ln„(/), as given in Eq. (|1.2p . Here, Eq. (|3.10l) . 
has been obtained by employing a physical mechanism 
based, exclusively, on the Lorentz transformations and 
without invoking any extra principle or assumption. 



As conclusions we recall the main results obtained in 
the present Letter: 

(i) We have shown that the Lorentz transformations 
are able to furnish a very transparent theoretical support 
to the function exp^(a;). This function results to be a 
light cone variable and therefore factorizes the energy- 
momentum Lorentz transformations. 

(ii) We have shown that the mathematical property 
exp^(a;) exp^(— cc) = 1 is enforced exclusively by the 
Galilei relativity principle. 



(iii) Eq. (|3.10l) . introduced heuristically in the lit- 
erature after 2002, whose solution defines the function 
In^ (x) , has been obtained here starting directly from the 
Lorentz transformations. 
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